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Abstract. We investigate the internal rotation profile of the GW Vir (PG1 159-type) 
star PG 0122+200 by employing an asteroseismological model that closely reproduces 
the observed periods of this star. We adopt a forward approach and two inversion meth- 
ods based on the rotational splitting of the pulsation frequencies to explore the proper- 
ties of the rotation of PG 0122+200. We found evidence for differential rotation in this 
star. 



1. Introduction 

In recent years, asteroseismology has become a powerful tool to unravel the internal 
structure of oscillating stars, in particular, pulsating white dwarfs and pre-white dwarfs 
(Winget & Kepler 2008, Althaus et al. 2010). The approach is suitable also to place 
constraints on the stellar rotation. In fact, stellar rotation removes the intrinsic mode 
degeneracy of a nonradial g-mode characterized by an harmonic degree £ and a radial 
order k. As a result, each pulsation frequency is split into multiplets of (2£ +1) fre- 
quencies specified by different values of the azimuthal index m, with m — 0, ±1, . . . , ±€ 
(Unno et al. 1989). If the rotation rate is slow compared with the pulsation frequencies, 
the frequency separation between each component of the multiplet is proportional to 
the rotation velocity of the star (Q). In this work, we perform a detailed asteroseismo- 
logical study aimed at placing constraints on the internal rotation of the GW Vir star 
PG 0122+200, using the best existing evolutionary and seismic model for this star. 



2. PG 0122+200 and the asteroseismological model 



PG 0122+200 has T eS = 80000 ± 4000 K and logg = 7.5 ± 0.5 (Dreizler & Heber 
1998). This pulsating PG1 159 star currently defines the locus of the low-luminosity red 
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edge of the GW Vir instability strip (Althaus et al. 2010). In this paper, we employ the 
high-quality observational data on PG 0122+200 gathered by Fu et al. (2007), consist- 
ing of 23 frequencies corresponding to modes with £ - I that include seven rotational 
triplets (m = -1,0, +1) and two isolated frequencies with (probably) m = — see 
Table 5 of Fu et al. (2007). The theoretical frequency splittings were assessed employ- 
ing the non-rotating asteroseismological model of PG 0122+200 derived by Corsico et 
al. (2007). This model reproduces the m - observed periods with an average of the 
period differences (theoretical vs. observed) of < 0.9 s. Employing an asteroseismo- 
logical model represents a substantial improvement over previous works of this kind — 
see, e.g., Charpinet et al. (2009), for the case of PG 1159-035. 

3. The forward approach 

This method has been described in detail in Charpinet et al. (2009). The theoretical 
frequency splittings are obtained varying the assumed rotation profile and then they are 
compared with the observed ones until a best global match is found. The theoretical 
splittings are computed using the expressions resulting from the perturbative theory to 
first order in Q, that assumes that the pulsating star rotates with a period much longer 
than any of its pulsation periods (Unno et al. 1989. The goodness of the match between 
theoretical (6vJ) and observed (<5v9) rotational splittings is described using a quality 
function defined as: 



Each difference is weighted with the inverse of the squared standard uncertainty (cr,) of 
the observed splittings, which is derived from the uncertainties in the frequencies given 
in Fu et al. (2007). 

We used very simple functional forms for £l(r), and in particular, linear differential 
rotation profiles, H(r) = (Q s - Q c ) r + f2 c , where O s and H c are the rotation rates 
at the stellar surface and center, respectively. We performed forward computations 
varying the parameters Q s and Q. c in the range 0-20 /iHz. We found that there exists 
a unique, well-localized best-fit solution at Q. c = 10.62 ± 1.8y«Hz and Q s ) = 4.41 ± 
1 . 1 /iHz. The quoted uncertainties are derived by means of an error analysis in which 
we assume that the uncertainties of the observed frequency splittings have a Gaussian 
distribution with a standard deviation of - 0.08 //Hz. The existence of this solution 
suggests that the central regions of PG 0122+200 could be rotating more than twice 
faster than the surface. These results disagree with those of Corsico & Althaus (2010), 
that indicate that the surface of PG 0122+200 could be spinning faster than the core. 
The discrepancy has its origin in the fact that, in that preliminary work, the fits of the 
rotational splittings were made without weighting the terms of the sum in the quality 
function, and thus, the impact of the different uncertainties of the observational data on 
the final result was neglected. 

Vauclair et al. (201 1) have found that seven oscillation frequencies of PG 0122+200 
are experiencing changes over time, with much larger amplitudes and shorter time 
scales than those expected by cooling. We have redone our forward computations by 
considering these frequency drifts. We have arrived at the conclusion that, even with 
the uncertainties produced by the changes of some frequencies of PG 0122+200 over 
time, rigid rotation can be discarded. 
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Figure 1. A test of our RLS method. The dashed straight line shows the input 
rotation profile (Q c = 1CT 4 //Hz, fi c = 20 //Hz), and the thin solid curve corresponds 
to the inverted profile for increasing values of the parameter A. 



4. Other inversion methods 

We have gone beyond the forward approach, and employed another different, indepen- 
dent method to explore the rotation profile of PG 0122+200. This is the Regularized 
Least Squares (RLS) fitting technique — see Kawaler et al. (1999) for details. We first 
tested the reliability of our RLS scheme by employing this technique on "synthetic" 
frequency splittings generated with the asteroseismological model of PG 0122+200 
through the forward approach. Specifically, we considered rotational splittings corre- 
sponding to consecutive £ = 1 g-modes with k - 1 , • • • , 40. The regularization matrix 
corresponded to the smoothing of the second derivative of O(r). In all of the cases ex- 
amined so far, the inversions were able to recover the input rotation profile we used to 
compute the synthetic splittings, provided that an adequate range of values of the reg- 
ularization parameter A was adopted. Fig. Q] illustrates the particular case of a linearly 
increasing rotation profile. Next, we applied the RLS method to infer the internal rota- 
tion profile of PG 0122+200. We have employed the seven observed (averaged) I = 1 
splittings. In Fig. [2] we show the inverted rotation profiles. For very small values of 
A, the inverted profiles exhibit strong variations that lack physical meaning. However, 
as the value of A is increased, the inverted solution gradually stabilizes. The resulting 
rotation profile (corresponding to A > 10~ 2 ) consists of an almost linearly decreasing 
rotation rate with O c ~ 10.75 //Hz and O s ~ 4.58 //Hz, in excellent agreement with the 
results of the forward approach. An analysis of errors similar to the one performed for 
the forward approach leads to the conclusion that even with the inclusion of reasonable 
uncertainties in the observed splittings, the rotation of PG 0122+200 is faster at the 
central regions than at the surface. Specifically, we found Q c = 10.75 ± 2.4//Hz and 
O s = 4.58 ±1.7 //Hz. 

Finally, we have also made rotational inversions onto a fixed functional basis. The 
method is called "function fitting" — for details, see Kawaler et al. (1999). Specifically, 
we explored linear functional forms for O(r). The optimal values we found are f2 c = 
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Figure 2. Inverted rotation profiles for PG 0122+200 corresponding to different 
values of the regularization parameter A. 



10.74 ± 2.9/iHz and Q s = 4.57 ± 1.8 //Hz, in excellent agreement with the RLS fits and 
also with the forward approach. 

5. Conclusions 

The three methods employed in this work to explore the internal rotation of PG 0122+200 
suggest that it is differential, with the central regions rotating more than twice faster 
than stellar surface. This constitutes the first asteroseismic evidence of differential ro- 
tation with depth for an evolved star. 
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